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SPECIFICATION 


TO ALL WHOM IT MAY CONCERN: 

Be it known that I, Mats A. Brenner, a citizen of Sweden, 
residing at 4275 Deerwood Lane, in the County of Hennepin and State of 
Minnesota have invented a new and useful APPARATUS FOR NAVIGATION 
SATELLITE SIGNAL QUALITY MONITORING, of which the following is a 
specification. 


APPARATUS FOR NAVIGATION SATELLITE 
SIGNAL QUALITY MONITORING 

Technical Field of the Invention 

The present invention relates generally to 
5 satellite based positioning systems such as the Global 
Positioning System (GPS) and more particularly to the 
monitoring of the quality of the signals transmitted by 
satellites in a satellite based positioning system. 

Background of the Invention and Prior Art 
10 A satellite based positioning system is used 

to determine a position of a receiver and typically 
includes satellite control facilities, a plurality of 
satellites, the receiver, and one or more local or 
regional ground stations. Each of the satellites 
15 transmits a signal that contains a code and certain 
prescribed information useful to the receiver in 
determining its position. The receiver synchronizes 
itself to the codes of at least four satellites and uses 
the information in the signals from these satellites in 
20 order to perform a triangulat ion like procedure so as to 
determine its coordinates and time offset with respect to 
a reference, such as the center of the Earth and the GPS 
standard time. 
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The receiver is not constrained to a specific 
location and, therefore, represents a variable position. 
Indeed, the purpose of the satellite based positioning 
system is to make it possible for the receiver to 

5 determine its position regardless of the location of the 
receiver. On the other hand, the local or regional 
ground station is in a fixed location and is used to 
monitor the signals transmitted by the satellites. The 
signals transmitted by the satellites can be adversely 

10 affected, for example, by atmospheric conditions which 
can lead to improper position determinations by the 
receiver. The ground station, therefore, notifies the 
receiver of any necessary signal corrections to allow the 
receiver to make more accurate position calculations. 

15 This arrangement is referred to as differential 

positioning . 

The ground station of the present invention 
also monitors the signals transmitted by the satellites 
in order to detect faults within the satellites. For 
20 GPS, these faults are specified by the FAA who imposes 
stringent requirements to protect users against 
positioning system signal faults. A set of test 
waveforms has been chosen by the FAA to represent at 
least some of the more egregious faults. These waveforms 
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are used for certification testing of the ground station 
equipment . 

The prior art determines faults by comparing 
conventional code tracking discriminators at different 
correlator spacings . As shown in Figure 1, a correlation 
curve is established by correlating the code received 
from a satellite with a suite of code references which 
are time shifted replicas of the code transmitted by that 
satellite. For example, seven correlation measurements 
may be calculated as shown in Figure 1. The in-phase 
measurement IP represents the amount of correlation 
between the received code and a reference code that has a 
zero time shift with respect to the received code (this 
measurement is referred to as punctual) . The in-phase 
measurement IE X represents the amount of correlation • 
between the received code and a reference code that has a 
first predetermined time shift so that it is early with 
respect to the received code. The in-phase measurement 
ILi represents the amount of correlation between the 
received code and a reference code that has a second 
predetermined time shift so that it is late with respect 
to the received code. Similarly, the in-phase 
measurement IE 2 is derived using a third predetermined 
time shift, the in-phase measurement IL 2 is derived using 
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a fourth predetermined time shift, the in-phase 
measurement IE 3 is derived using a fifth predetermined 
time shift, and the in-phase measurement IL 3 is derived 
using a sixth predetermined time shift. The magnitude of 
the first predetermined time shift may be equal to the 
magnitude of the second predetermined time shift, the 
magnitude of the third predetermined time shift may be 
equal to the magnitude of the fourth predetermined time 
shift, and the magnitude of the fifth predetermined time 
shift may be equal to the magnitude of the sixth 
predetermined time shift. It is assumed that all 
measurements are normalized such that the measured 
correlation is a function of the time shifts only and not 
the absolute power of the received satellite signal. 

First, second, and third discriminators are 
then formed according to the following equations: 


d x = {IL X - IE X ) IP 
d 2 = (IL 2 - IE 2 ) IP 
d 3 - (IL 3 - IE 3 ) IP 
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These discriminators are thereafter compared to each 
other through the formation of quantities d 12 , d 13 , 
d 2 3 according to the following equations: 

d lt2 = K - eg 


d i, 3 = K " d 3 | 

d 2,3 = l d 2 " d al 

The quantities d lf2 , d 1#3 , and d 2>3 are compared to 
corresponding thresholds D 1#2 , D lf3 , and D 2>3 such that, if 
the first discriminator d lf2 exceeds D 1<2 , if the second 
discriminator d lj3 exceeds D lf3 , or if the third 
discriminator d 2#3 exceeds D 2<3 , a fault is assumed to 
exist. During normal operation of the global positioning 
system, this test is performed on the signals received 
from each of the satellites. During certification, a 
test is to be performed using each of the test waveforms 
chosen by the FAA in order to prove that fault detection 
occurs . 

At least one of the problems with this method 
is that it is requires six correlators in order to 
determine the three quantities d lf2 , d x , 3 , and d 2/3 which is 


Attorney Docket 
A13-26546 


too much hardware for the amount of useful data being 
provided . 

It is also known for ground stations to 
determine faults by scanning the whole correlation peak 
5 (i.e., the portion of the correlation curve around the 
punctual in-phase measurement IP) in order to determine 
whether the peak varies from some prescribed norm by a 
predetermined amount. However, this fault detection 
arrangement requires a substantial amount of computing 
10 power and it lacks accuracy. 

A third method in the prior art uses the 
following ratios between the measurements IE3, IE2 , IE1, 
IL1, IL2 , and IL3 : 


IE3 


15 


■E3.E2 j E2 

IE3 

-E3,E1 j E1 

IE3 

C E3,L1 j Ll 
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Each of these ratios is compared to a corresponding 
predetermined value . 

The present invention is directed to an 
arrangement which overcomes one or more problems of the 
5 prior art . 

Summary of the Invention 

In accordance with one aspect of the present 
invention, an apparatus for the detection of positioning 
system satellite signal faults comprises a correlator and 

10 a fault detector. The correlator determines a plurality 
of correlation measurements at points along a correlation 
curve, and each correlation measurement is based upon a 
correlation between a received satellite signal and a 
reference. The fault detector determines differences 

15 between the correlation measurements along the 
correlation curve and detects a fault from the 

differences . 

In accordance with another aspect of the 
present invention, a method of detecting faults affecting 
20 a signal transmitted by a positioning system satellite 
comprises: correlating the transmitted signal with a 
first reference in order to determine a first correlation 
measurement at a first point along a correlation curve; 
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correlating the transmitted signal with a second 
reference in order to determine a second correlation 
measurement at a second point along the correlation 
curve; correlating the transmitted signal with a third 
reference in order to determine a third correlation 
measurement at a third point along the correlation curve; 
determining a first difference from the first and second 
correlation measurements; determining a second 
difference from the second and third correlation 
measurements; directly comparing the first difference to 
a first threshold; directly comparing the second 
difference to a second threshold; and, detecting a fault 
in the satellite based upon the comparisons of the first 
and second differences to the first and second 
thresholds . 

In accordance with still another aspect of the 
present invention, a method of detecting faults affecting 
a signal transmitted by a positioning system satellite 
comprises: correlating the transmitted signal with 
references in order to determine a plurality of 
correlation measurements at corresponding points along a 
correlation curve; determining a single value from n 
pairs of the correlation measurements, wherein n > 2; 
comparing the single value to a threshold; and. 
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detecting a fault in the satellite based upon the 
comparison . 

Brief Description of th e Drawings 

These and other features and advantages will 
5 become more apparent from a detailed consideration of the 
invention when taken in conjunction with the drawings in 
which : 

Figure 1 is a waveform showing a correlation 
diagram useful in explaining prior art fault detection as 
10 implemented in a ground station in a global positioning 
system; 

Figure 2 is a schematic diagram of a portion 
of a ground station receiver pertinent to the present 

invent i on ; and , 
15 Figure 3 is a waveform showing a correlation 

diagram useful in explaining fault detection as 
implemented by a ground station in a global positioning 
system in accordance with the present invention. 

Detailed Description 

20 A portion of a ground station 10 pertinent to 

the present invention is shown in Figure 2. The ground 

station has correlators 12-E m/ . . . , 12-E 3 , 12-E 2/ 12-E lt 
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12-P, 12-L 1# 12 -L 2 , 12-L3, . . 12-L n/ where n + m is 

greater than two, and where n is the number of late 
correlation measurements and m is the number of early 
correlation measurements to be used in determining a 
fault. The correlator 12-P correlates the usual code in 
the received signal with a reference 14 -P to produce a 
punctual correlation output IP, the correlator 12-L lf 
correlates the code in the received signal with a 
reference 14 -L x to produce a late correlation output IL 1# 
the correlator 12 -L 2 correlates the code in the received 
signal with a reference 14 -L 2 to produce a late 
correlation output IL 2 , the correlator 12 -L 3 correlates 
the code in the received signal with a reference 14 -L 3 to 
produce a late correlation output IL 3 , . . and the 

correlator 12 -L n correlates the code in the received 
signal with a reference 14 -L n to produce a late 
correlation output IL n - 

In addition, a correlator 12 -E x correlates the 
code in the received signal with a reference 14 -Ex to 
produce an early correlation output IE 1# a correlator 12- 
E 2 correlates the code in the received signal with a 
reference 14 -E 2 to produce an early correlation output 
IE 2 , a correlator 12 -E 3 correlates the code in the 
received signal with a reference 14 -E 3 to produce an 
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early correlation output IE 3 , . . ., and a correlator 12- 
E m correlates the code in the received signal with a 
reference 14 -E m to produce an early correlation output 

5 The ground \stat ion 10 has a processor 16 which 

uses the punctual and Ate correlation outputs IP, IL 1# 
IL 2 , IL 3 , . . ., IL n as disclosed hereinafter in order to 
determinei whether a faultt exists. Alternatively or 
additionally, the processor\l6 can use the early 


10 


correlation outputs IE 1# IE 2 ,\lE 3 , . . . , IE m as disclosed 
hereinafter in order to detert^Lned whether a fault 
exists . 

In order to generate the punctual correlation 
output IP, the processor 16 shifts the reference 14-P, 
15 which may be a replica of the code contained in the 
received signal, until an optimum correlation is 
obtained. The processor 16 then controls the reference 
14-L-l so that the reference 14-L X is a replica of the 
reference 14-P and so that the reference 14 -L x is time 
shifted with respect to the reference 14-P by a first 
predetermined amount of time. Accordingly, the 
correlator 12 -L x produces the late correlation output IL X 
The processor 16 also controls the reference 14 -L 2 so 
that the reference 14 -L 2 is a replica of the reference 


20 
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14 -P and so that the reference 14-L 2 is time shifted with 
respect to the reference 14 -P by a second predetermined 
amount of time, where the second predetermined amount of 
time is greater than the first predetermined amount of 
5 time. Accordingly, the correlator 12 -L 2 produces the 

late correlation output IL 2 . Similarly, the processor 16 
controls the reference 14 -L 3 so that the reference 14 -L 3 
is a replica of the reference 14-P and so that the 
reference 14 -L 3 is time shifted with respect to the 
10 reference 14-P by a third predetermined amount of time, 
where the third predetermined amount of time is greater 
than the first and second predetermined amounts of time. 
Accordingly, the correlator 12 -L 3 produces the late 
correlation output' IL 3 . The remaining late correlation 
15 outputs up to IL n are generated in a like manner. The 

first, second, third, etc. predetermined amounts of time 
are all chosen so that the late correlation outputs IL X 
through IL n are all on the downward or late slope of the 
correlation curve as shown in Figure 3 . 
20 Additionally or alternatively, the correlators 

12-Ei, 12-E 2 , 12-E 3 , . . ., 12-E m may be positioned so as 
to generate the early correlation outputs IE lf IE 2 , IE 3 , . 

IE m . Also, quadrature phase correlation outputs QE m , 
. . . , QEi, QP, QL X , . . . , QL n may be generated by 
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correlating the code in the received signal to a time 
shifted quadrature form of the reference 14-P. In 
accordance with this latter alternative, each measurement 
used to generate a fault indication may be formed as an 
5 RMS (Root Mean Square) value of the corresponding in 
phase and quadrature phase measurements. 

The set IE m , . . ., IE 3 , IE 2 , IE X , IP, IL 1# 
IL 2 , IL 3 , . - ., IL n may be denoted as I m , . . ., I- 3 , z -2> 
I.i, I 0 » 1 i> T 2< Is/ • • ■> x n and the following 
10 corresponding set of RMS values 


15 


20 


may be denoted as R m , . . ., R. 3 ' R -2' R -i» R o> R i' R 2' R 3' 
• • ■ 1 R n - 

If early as well as late correlation outputs 
are to be used for fault detection, the processor 16 
processes the early correlation outputs IE m through IE lf 
the punctual correlation output IP, and/or the late 
correlation outputs IL3. through IL n so as to derive one or 
more measured differences d i(j . These measured 
differences d it j are generated in accordance with the 
following equations: 
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d. . = I, - I, <!> 

if 3 i J 


or 


d. .= !?,- (2) 


5 where i = -m, . . . , n and j = -m, . - . , n, and where 
the negative sign indicates measurements on the early 
slope and the positive sign indicates measurements on the 
late slope of the correlation curve. 

At this point, it is possible to subtract the 

10 expected difference from all or a subset of these 

measured differences d irj and to compare the resulting 
difference deviations to corresponding thresholds in 
order to determine the existence of a fault. For 
example, assuming that all of these difference deviations 

15 are used, then these difference deviations may be 

compared to corresponding thresholds in accordance with 
the following equation: 


(3) 
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where Ed i>;j is the difference that is expected for each 
corresponding measured difference d i(j when there is no 
fault . 

In some cases, the measured differences d i(j 
5 may be affected by thermal and multipath noise which 

could lead to false detection of faults, depending upon 
the sensitivity of the fault detection apparatus, i.e., 
the magnitudes of the thresholds D i;j . Accordingly, in 
these cases, a fault could be detected when no fault is 
10 in fact present, or a fault which is present might not be 

detected at all. 

The thermal noise content in d i :j can be 
determined as a function of the delay h ifj between the 
reference codes 14-E m , . . ., 14-E 3 , 14-E 2 , 14-E^ 14-P, 

15 14 -L x , 14-L 2 , 14-L 3 14-L n . The delay h^ is the 

delay between the two references that are correlated with 
the received signal to produce Ii and Ij . Typically, h if j 
= 0.025 to 0.05 chip, but may vary from this range. The 
thermal noise thl in d ifj depends on the signal to noise 

20 ratio and the standard deviation (1-sigma) of thl and is 
given by the following equation: 


h. .B 

a (i i) = 293. 1,J ( 4 ) 
°thi^ 1 '^ ) s/No 
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where B is the two-sided bandwidth of the noise. In 
addition, there is another contribution, th2 , to the 
thermal noise due to the variation of the punctual 
reference (i.e., the reference 14-P) . Accordingly, the 
total thermal noise is th = thl + th2 . The multipath 
noise mp depends on the antenna gain pattern and its 
overbounding 1-sigma a^d^) ( "e ) is expressed as a 
function of satellite elevation *e . The statistical 
properties of th and mp are usually identified at 
installation of the ground station and the statistical 
information is parameterized and are thereafter stored in 
memory . 

One way to minimize any adverse effects of 
thermal and multipath noise is to make a plurality of 
measurements for each of the measured differences d i#j 
that are used in the detection of faults. Then, the 
measurements for each of the measured differences d i#j may 
be averaged or filtered. Because the thermal noise and 
some of the multipath noise are not particularly 
correlated from one measurement to the next, averaging 
will tend to reduce the effects of thermal and multipath 
noise . 

As an example, let it be assumed that the 
punctual correlation output IP and the late correlation 
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outputs IL X and IL 2 are used to detect faults. 
Accordingly, the following measured differences are 
determined: d 0>1 = IP - 11^; d 0<2 = IP - IL 2 ; and, d 1/2 = 
IL X - IL 2 . In order to reduce the effects of thermal and 
multipath noise, however, plural calculations of the 
measured difference d 01 are made based upon plural 
correlation measurements resulting in plural punctual 
correlation outputs IP and plural late correlation 
outputs ILi- All such calculations of the measured 
difference d 0<1 are then averaged. Similarly, plural 
calculations of the measured difference d 0<2 are made 
based upon the plural correlation measurements resulting 
in plural punctual correlation outputs IP and plural late 
correlation outputs IL 2 . As before, all such 
calculations of the measured difference d 0<2 are averaged. 
Likewise, plural calculations of the measured difference 
d 1>2 are made based upon the plural correlation 
measurements resulting in the plural late correlation 
outputs IL X and plural late correlation outputs IL 2 . 
Again, all such calculations of the measured difference 
d x 2 are averaged. These averages may then be compared to 
their corresponding thresholds D 0il , D 0<2 , and D 12 in order 
to determine the existence of a fault. 
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Another way to reduce the effect of thermal 
and multipath noise is to suitably filter the measured 
differences d i;j , or the punctual correlation output IP, 
the late correlation outputs IL X through IL n , and the 
5 early correlation outputs IE X through IE m , such as with a 

low pass filter. 

Still another way to reduce the effect of 
thermal and multipath noise is by implementing the 
following procedure. In describing this procedure, it is 
10 useful to define a covariance matrix P in accordance with 
the following equation: 

P = E[(d - m) (d - m) T ] < 5 ) 

where the underlines indicate vectors, where E [A] is the 
statistical expectation of A, where the vector m is the 
15 mean value of the vector d, and where the vector d is 
determined in accordance with the following equation: 


d T = {d ir d 2 ,d 3 ,d i , . . . ,d N ) 


(6) 


where d k = I k - 1^ - Ed k for k = -m n-1 or where 

dk = Rk _ Rk _ x - Ed k for k = -m n-1 assuming N + 1 

20 correlation measurements such as I m , . . ., I_ 3 , 1-2/ 
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I 0 , l 1# I 2 , I 3 I n . An upper triangular matrix U 

and a diagonal matrix D are determined according to the 
following equation: 

P = UDU T (7) 

where P is the covariance matrix given by equation (6) . 
With the covariance matrix P known from equation (6) , the 
upper triangular matrix U and the diagonal matrix D can 
be determined, for example, by using Cholesky 
factorization. Thus, the following relationship may be 
defined in accordance with the following equation: 

d = CT 1 (d - m) (8) 

where d is a vector representing the decorrelated 
deviations generating the vector d. Equation (9) can be 
re-written according to the following equation: 

d = U d + m O) 


Then, combining equations \6) and (10) produces the 
following equation: \ 
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P= E [Ud [Ud) T ] = UE [d (d) r ] U T 


(10) 


By comparing equatiori\s (6) and (11), it can be seen that 
D is given by following equation: 


D = E [d (d) T ] 


(11) 


5 and that D, as defined above, is a diagonal matrix having 
the following format: 


D = 


a 0 0 ... 0 
0 o 0 ... 0 
0 0 o ... 0 


10 0 0 


(12) 


10 


Variances o 2 ± are then determined from the diagonal matrix 


D. 


As can be seen from the above equations, the 


deviations in the vector d where i varies from 1 to N 


are uncorrelated and have the variances o 


~2 


20 
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The final x 2 value for determining a fault is 
obtained according to the following equation: 


n d 2 


d[x 2] = j .2* (13) 


A normalization to o\ 1 as required in the definition of 
5 x 2 will be performed i\ equation (14). The value d[x 2 ] is 
a single value which ha sV reduced thermal and multipath 
noise, which represents information regarding a plurality 
of correlation measurement! , and which may be compared to 
a threshold D in order to determine the existence of a 
10 fault. 

Certain modifications of the present invention 
have been discussed above. Other modifications will 
occur to those practicing in the art of the present 

invention. For example, as described above, the x 
15 distribution is based on the assumption that all involved 
distributions are Gaussian. The distributions of d k may 
deviate from this assumption and appropriate corrections 
to the formulas given here may be necessary. 

Moreover, the present invention has been 
described above in connection with the detection of 
satellite signal faults such as those specified by the 


10 
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FAA. These faults result in signal distortions 
detectable by use of the present invention. The present 
invention as embodied by the following claims can also be 
used to detect other signal distortions such as those 
arising from multipath and satellite code cross 
correlation effects. 

Accordingly, the description of the present 
invention is to be construed as illustrative only and is 
for the purpose of teaching those skilled in the art the 
best mode of carrying out the invention. The details may 
be varied substantially without departing from the spirit 
of the invention, and the exclusive use of all 
modifications which are within the scope of the appended 
claims is reserved. 


